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Ion conducting {CaO) x (S1O2) \- x glasses and melts show a threshold be- 
haviour in dc conductivity near x = Xt = 0.50, with conductivities increasing 
linearly at x > x%. We show that the behaviour can be traced to a rigid 
(x < 0.50) to floppy (x > 0.50) elastic phase transition near x = x%. In 
the floppy phase, conductivity enhancement is traced to increased mobility or 
diffusion of Ca 2+ carriers as the modified network elastically softens. 
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Electrical and thermal properties are directly related to the structure of silicate melts 
and investigations into their structure have received therefore considerable attention in earth 
science. It is indeed of fundamental importance to understand how magmas which contain 
alkali and alkaline earth silicates behave under different physical circumstances. Alkali 
silicate glasses have been extensively studied within this context and numerous articles 
have focused on the variation of structure [1], conductivity [2], glass transition [3], etc. 
with respect to the alkali concentration. On the other hand, a small number of studies on 
alkaline earth silicate glasses have been reported to date. Some of them, however, have 
focused either on the structure of glasses containing barium and calcium silicates [4], or on 
the thermodynamics and miscibility of the corresponding melts [5] and density or specific 
volume [6] . It is worth mentioning that most of these studies have been restricted around the 
50% concentration of alkaline earth oxide where glass-forming tendency (GFT) is optimized 
[7]. Finally, even though most of the fast ionic conductors (FIC) [8,9] having potential 
applications in electrochemical devices such as solid state batteries are alkali (or silver) 
silicates or thiosilicates, the conductivity of alkaline earth silicates is non-negligible [about 
2.10~ 7 n^.cm^ 1 at 550°C] and the role of the migrating calcium cations remains to be 
completely understood. Our work attempts to address for the first time this basic issue. 

In this Letter, we report on (CaO) x (Si02)i- x glasses and the common physical origin 
for the behavior in both the GFT and the conductivity data that follow directly from the 
Phillips-Thorpe constraint theory [10]. Calcium silicate melts and glasses near x = 0.5 are 
optimally constrained and display percolative conduction in both melts and glasses. This 
is a new feature that has never been observed in the corresponding alkali silicates, which 
links ionic conductivity to glass elasticity and structure. Furthermore, we show that the 
sudden increase of the conductivity at high calcium concentration has to be a consequence 
of the calcium cation mobility percolation arising from the floppy to rigid transition [11]. The 
major consequence is that here the carrier concentration does not dominate the conductivity, 
as usually reported. 

The samples were prepared by mixing pre-dried Si0 2 (99.99 %) and CaC0 3 (99.95 %) 
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powders in the correct proportions. For each sample, the mixture was melted in a platinum 
crucible at 1650°C for two hours, and quenched by plunging the bottom of the crucible in 
cold water. The glass transition temperatures were determined with a differential scanning 
calorimeter Setaram DSC-1600 at a heating rate 10°C/min. These values are slightly larger 
than those reported by Shelby [12] using the dilatometric technique but they exhibit the 
same global trend (fig. 1), i.e. a plateau up to a Calcium concentration in the range [0.45- 
0.48] and a more or less linear increase for higher concentrations. In the solid state, the 
complex electrical conductivity was measured with an impedance spectrometer 4194A in 
the temperature range 500°C — 950°C and in the frequency range 100 Hz-10 MHz. The 
samples were 0.5 cm 2 glass pellets with a thickness about 1 mm. Platinum was evaporated 
as electrodes on both faces of the pellet. The temperature was measured with a Pt/Pt- 
10%Rh thermocouple located at about 1 mm from the sample [13]. We restricted our study 
in the molten state to only three compositions. 

The dc conductivity a displays almost the same trend as T g (x), with a very small value 
(about 10~ 7 Q^.cm^ 1 ) for x < 0.47 and a sudden increase for larger compositions. The 
almost constant value of a in the low calcium region and the sudden increase for x > 0.47 has 
been also obtained in simulations and experiments for the molten state (fig. 2, insert [14]), 
suggesting that molten and glassy state behave similarly. Furthermore, we have observed 
the same kind of trend in our molten samples (see fig. 3, below) which are almost at the 
same temperatures than those of Ref. [14]. Figure 3 shows the Arrhenius plots [oT = 
do exp[— E^/hsT]] of the conductivity a for the glassy and the molten state. From the 
latter, it is also obvious that the conductivity data of three compositions in the x < 0.50 
range map onto each other and have the same global trend in the glass. 

These results can be quantitatively understood using the tools of Lagrangian bonding 
constraints as introduced by Phillips [10]. A network constrained by bond-stretching and 
bond-bending forces sits indeed at a mechanically critical point when the number of con- 
straints per atom n c equals the network dimensionality. At this point, glass optimum is 
attained and GFT is enhanced [15]. These ideas can be cast in terms of percolation theory 
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by evaluating the number of zero frequency modes (floppy modes F/N) from a dynamical 
matrix [11], which vanish at the floppy to rigid transition. The very accurate agreement of 
experiments in chalcogenides with the predictions of the theory is quite remarquable [16,17] 
but further investigations devoted to alkali oxide glasses have been performed only recently 
[18,19]. 

We consider the CaO — S1O2 system as a network of N atoms composed of n r atoms that 
are r— fold coordinated. Enumeration of mechanical constraints gives for bond-stretching 
forces r/2 and for bond-bending forces (2r — 3) constraints. The average number of floppy 
modes per atom F/N in a three-dimensional network is given by [11]: 

1 _ )V 

F/N = n d -n c = 3--Y,n r (--3) (0.1) 

iV r>2 Z 

which applied to the present (CaO) x (Si02)i- x case, yields: 

11 — 7x , 
F/N = 3- (0.2) 

provided that Si is four-folded and calcium and oxygen are two-fold coordinated. The 
coordination number used in equ. (0.2) for Ca desserves some comments since computer 
simulations [14] and EXAFS studies [20] have suggested that the number of nearest neigh- 
bors of the calcium atom may be about 6. However, Debye- Waller factor in these studies 
on combined systems using either BaO or CaO as network modifier are slightly different, 
suggesting that the oxygen neighbors of Ca were not all equivalent [20] thus promoting the 
chemical hypothesis for the coordination number of Ca (r = 2). On the other hand, in con- 
straint theory one deals only with the bonds that provide a strong mechanical constraint, 
which means that the covalent bonds are the most qualified. To be more specific, one can 
use Pauling's definition to determine the fractional ionic character / in the Ca-0 bond [22], 
which yields here the value / = 0.774 and a corresponding covalency factor of 0.226. With 
a coordination of 6 for Ca, the corresponding covalent coordination is 1.36, somewhat lower 
than 2. Furthermore, mechanical constraints provided by alkali atoms connected to more 
than one oxygen have been suggested to be only resonating constraints [19], which do not 



participate in the global counting. Finally, alkali silicate and tellurate glasses are illustrative 
precedents [18]. In these systems, the floppy to rigid transition occurs at the expected criti- 
cal concentration [21] only if the involved atoms (Na) have a covalent coordination number 
of 1 (and the "Pauling" analysis [23] yields the covalent coordination 0.832 — 1.04), even 
though the number of nearest neighbors [24] of the sodium atom is believed to be 5. 

Our interpretation of the dc conductivity is as follows. At x = x c = 0.5, the number of 
floppy modes F/N vanishes and the network undergoes a rigid to floppy transition, following 
equ. (0.2). For x < 0.5 the system is stressed rigid, i.e. there are more constraints than 
degrees of freedom. As a consequence, the mobility of the calcium cation is very weak 
because the cation has to overcome a strong mechanical deformation energy to move from 
one anionic site to another. In an ideal floppy glass at x > 0.5 where only bond-bending and 
bond-stretching forces are considered [10] this deformation energy is zero and related elastic 
constants (en, C44) vanish [25]. Therefore, percolation of Soppiness equals percolation of Ca 
mobility. As a result, one has a substantial increase of the mobility hence of the conductivity. 
In the strong electrolyte Anderson-Stuart model [26], the activation energy for conductivity 
E A = E c + E m depends on the Coulombic term E c (which controls the carrier rate) and 
a strain term E m related to the mobility that can be thought as the energy required to 
enlarge the radius of an anionic site of length l(x) by an amount Sr, the latter quantity 
being roughly equal to the cation radius: 

E m (x) = ^irc 44 (x)l(x)(5r) 2 (0.3) 

The energy term E m (x) depends on the shear modulus 044(0;) which in turn depends on 
the Ca concentration [29]. In a floppy model network, the shear modulus 044(0;) is zero [25] 
leading to a zero value for the strain energy E m (x) and a maximum in mobility. 

Additional support for our interpretation derives from two other observations: (1) In the 
molten (CaO) x — {SiO^i-x-, we still observe an increase in conductivity with x but with 
values which are substantially higher compared to the glass (fig. 3) and which can be fitted 
with a Vogel-Fulcher-Tamman (VFT) law: 
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oT = Aexp[-A/(T - T )] (0.4) 

arising from the variation in temperature of the calcium diffusion constant (via the Nernst- 
Einstein equation) [26]. From the fits (solid lines in insert of fig. 3) it appears that the 
pseudo-activation energy A is constant for the two values in the floppy region (A = 2733.0 
and 2740. 7K for respectively x = 0.44 and x = 0.5) while it increases for the composition 
x = 0.53 (A = 3008 K). Also, the VFT temperature T at which the diffusion constant 
(and the underlying relaxation time towards thermal equilibrium) diverges is dynamically 
inaccessible during glass transition, but it is accepted that its value is close to the ideal 
Kauzmann [27] glass transition temperature T K . An ideal glass would therefore have its T g 
close to T K and T , a situation which is met in the present system for the Phillips' optimal 
glass composition x = 0.5 (T g /T = 1.0002 compared to T g /T = 1.0218 for x = 0.44 and 
T g /T = 1.0621 for x = 0.53). 

In fig. 4 are represented the Arrhenius activation energy for conductivity and the 
preexponential factor cr as a function of the mean coordination number of the network [11] 
which is f = (8 — Ax)/ (3 — x). The former displays a minimum at f = 2.4 which corresponds 
to x = 0.5. 

This in turn can be compared to the non-reversing heat flow extracted from complex 
calorimetric measurements at the glass transition temperature for various chalcogenide sys- 
tems [30]. Here, this relaxing part of the total heat flow exhibits always a minimum at the 
floppy to rigid transition (close to the mean- field [10] value f = 2.4). The same global trend 
is observed in numerous chalcogenides [30], and our calcium silicate glasses exhibit also a 
minimum in at the mean coordination number f = 2.4. From the chalcogenide exam- 
ple, we suggest that when x = 0.5 structural relaxation of the glass network proceeds with 
minimal enthalpic changes because the minuscule network stress is uniformly spread when 
the network is mechanically critical (i.e. F/N=0). Furthermore, we note that the variation 
of a is here not driven by the compositional trends of ctq. It appears to be rather a balance 
between the variation of Ea and a . (2) Rigidity percolation thresholds in Micro- Raman 
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measurements performed in backscattering geometry have been reported in Ref. [16]. There- 
fore, we have studied in the present work as a function of the Ca concentration the frequency 
v and the linewidth T of the stressed rigid S1O4/2 {Q A unit, where the superscript denotes 
the number of network bridging oxygens)). The corresponding line {A\ stretching mode) 
exhibits a change of regime (fig. 5) for the line frequency at the concentration x = 0.50, 
consistently with chalcogen analogs [16]. We concentrate here only on this (stressed rigid) 
line and will report separately the complete Raman analysis elsewhere [31]. On the other 
hand, the evolution of the linewidth T with x permits to follow the local environnement 
of the Q 4 unit. For x < 0.50, T remains constant, related to the absence of change in the 
coupling of this unit with the rest of the network. It is the coupling which makes possible 
the presence of rigid regions (through isostatic Q 4 — Q 3 and stressed Q 4 — Q A bondings) 
although the number of floppy Q 2 and Q l units is steadily increasing. Above the critical 
concentration x = 0.5, the sharp drop of T clearly shows decoupling of the Q A unit with re- 
spect to the network (fig. 5), signifying decoupling of stressed rigid regions thus percolation 
of Soppiness. 

In FICs, it is often believed that it is the carrier concentration that dominates the 
conductivity. In the strong electrolyte model [26], the cation electostatic Coulombic energy 
barrier has to be overcome to ensure conduction. On the other hand, in the weak electolyte 
model a dissociation energy is needed to create the mobile carrier [28]. These two pictures 
remain of course valid as long as the sizes of the cations is weak compared to the interstices 
of the glass network. 

Our conclusion brings us back to the analogy with alkali oxide FIC's and the popular 
conductivity channel picture [24]. In these glasses, the rigid to floppy transition occurs at 
the concentration x c = 0.20 which is very close to the reported threshold concentration 
separating intrachannel cation hopping (involving a weak mechanical deformation of the 
network, since the motion occurs only in macroscopic holes of the network) from network 
hopping (strong mechanical deformation only possible in a floppy network). However, in 
alkali silicates, no typical behavior emerges in compositional trends of the conductivity be- 
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cause of the growing contribution of the free carrier rate [24] . This has to be put in contrast 
with our present study on calcium silicates where the conductivity is almost constant up 
to a critical value x c ~ 0.48 beyond which a steadily increases. Therefore, the free carrier 
concentration cannot be considered as an increasing function of alkaline earth composition. 
We believe that conductivity in calcium silicates is driven by the carrier mobility in the net- 
work and percolates at the rigid to floppy transition. On this system, further investigations 
and measurements of physical quantities displaying usually a threshold behavior [33] at the 
rigidity transition will be achieved in close future. 

LPTL is Unite Mixte de Recherche n. 7600. It is a pleasure to acknowledge the help of 
Boris Robert during the course of this work. 
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FIG. 1. Compositional trends in glass transition temperature T g (x) in (CaO) x (Si02)i- x from 



DSC and electrical measurements, compared to previous studies[12] 
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FIG. 2. dc conductivity a in {CaO) x {Si02)\- x glasses with Ca composition x for different 
temperatures. The insert shows the conductivity in the melt from simulation and previous exper- 
iments^] at around 2000 K. 
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FIG. 3. Arrhenius plots of the ionic conductivity for (CaO) x (Si02)i- x glasses in the solid and 
the molten state. The fit for the estimation of Ea and ao have been performed at low temper- 
ature in the Arrhenius regime. Note the small deviation occuring close to the glass transition 
temperature. The insert shows oT in the liquid as a function of temperature with corresponding 
Vogel-Fulcher-Tamman (VFT) fits. 
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FIG. 4. Arrhenius parameters (activation energy Ea (solid circles) and preexponential factor cfq 
(open circles)) as a function of the mean coordination number f in the glass phase. 
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FIG. 5. A\ mode frequency v (open circles) and linewidth T (filled circles) of the Q 4 line as a 
function of Ca concentration. The insert shows part of the Raman spectra with the line of interest 
for x = 0.46 and the related Gaussian fits at 850 cvrT x [Q° unit] and 950 cm~ l [Q 2 unit] from [32]. 
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